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ABSTRACT 

We report the discovery of two new Galactic candidate luminous blue variable (cLBV) 
stars via detection of circular shells (typical of known confirmed and cLBVs) and 
follow-up spectroscopy of their central stars. The shells were detected at 22 /im in 
the archival data of the Mid-Infrared All Sky Survey carried out with the Wide-field 
Infrared Survey Explorer (WISE). Follow-up optical spectroscopy of the central stars 
of the shells conducted with the renewed Southern African Large Telescope (SALT) 
showed that their spectra are very similar to those of the well-known LBVs P Cygni 
and AG Car, and the recently discovered cLBV MN112, which implies the LBV clas- 
sification for these stars as well. The LBV classification of both stars is supported by 
detection of their significant photometric variability: one of them brightened in the 
R- and /-bands by 0.68 ± 0.10 mag and 0.61 ± 0.04 mag, respectively, during the last 
13-18 years, while the second one (known as Hen 3-1383) varies its B, V, R, I and K s 
brightnesses by ~ 0.5 — 0.9 mag on time-scales from 10 days to decades. We also found 
significant changes in the spectrum of Hen 3-1383 on a timescalc of ~ 3 months, which 
provides additional support for the LBV classification of this star. Further spectropho- 
tometric monitoring of both stars is required to firmly prove their LBV status. We 
discuss a connection between the location of massive stars in the field and their fast 
rotation, and suggest that the LBV activity of the newly discovered cLBVs might be 
directly related to their possible runaway status. 

Key words: line: identification - circumstellar matter - stars: emission-line, Be - 
stars: individual: Hen 3-1383 
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1 INTRODUCTION 



Massive stars go through a sequence of transitional evolu- 
tionary stages whose duration and order are still not well 
established (Chiosi & Maeder 1986; Langer et al. 1994; 
Stothers & Chin 1996; Maeder & Meynet 2011). Of these 
stages, perhaps the most important in the evolutionary sense 
and interesting in observational manifestations is the Lumi- 
nous Blue Variable (LBVs) stage (Conti 1984; Wolf 1992; 
Humphreys & Davidson 1994; van Genderen 2001; Vink 
2009). During this stage a massive star could exhibit ir- 
regular, strong spectroscopic and photometric variability on 
timescales from years to decades or longer. This variabil- 
ity manifests itself in changes of the stellar spectrum from 
that of late O/early B-type supergiants to that of A/F-type 
ones (e.g. Stahl & Wolf 1982; Stahl et al. 2001) and in the 
brightness increase by several magnitudes during the cooler 
(outburst) state. Major eruptions of LBV stars can masquer- 
ade as peculiar supernovae (Goodrich et al. 1989; Filippenko 
et al. 1995; Van Dyk et al. 2002; Smith et al. 2011), while 
some LBVs could perhaps be the direct progenitors of core- 
collapse supernovae (Kotak & Vink 2006; Smith et al. 2007; 
Gal- Yam & Leonard 2009). 

The spectrophotometric variability of LBVs is often ac- 
companied by drastic changes in the stellar mass-loss rate 
so that essentially all LBVs are surrounded by compact 
(~ 0.1 — 1 pc in diameter) shells (visible in the infrared, 
optical and/or radio wavebands) with a wide range of mor- 
phologies, ranging from circular to bipolar and triple-ring 
forms (e.g. Nota et al. 1995; Weis 2001; Smith 2007; Gvara- 
madze, Kniazev & Fabrika 2010c). The origin of the LBV 
shells can be attributed either to the interaction of stellar 
winds during the subsequent phases of evolution of mas- 
sive stars (e.g. Robberto et al. 1993; Garcia-Segura, Langer, 
Mac Low 1996) or to instantaneous ejections (e.g. Smith & 
Owocki 2006) caused by some catastrophic events (triggered 
either by internal or external factors). 

The LBV phenomenon is still poorly understood, which 
is mostly due to the scarcity of known LBV stars (caused by 
the short duration, a few 10 3 — 10 4 years, of the LBV phase). 
Several years ago there were only 12 confirmed Galactic 
members of the class and 24 candidate LBVs (cLBVs) 
(Clark, Larionov & Arkharov 2005; Vink 2009). The discov- 
ery of additional Galactic LBVs would, therefore, greatly ad- 
vance the field and help in understanding the evolutionary 
status and connection of LBVs to other massive transient 
stars, the diversity of the environments in which they form, 
and their distribution in the Galaxy. One can also expect 
that the increase of the number of known LBVs would allow 
us to understand the driving mechanism(s) of the LBV phe- 
nomenon and to answer the question: whether a deficiency 
of LBVs with luminosities between log(L / X©) ~ 5.6 and 5.8 
is real or is a result of the small-number statistics (Smith, 
Vink & de Koter 2004)? 

A detection of LBV-like shells can be considered as 
an indication that their associated stars are massive and 
evolved (e.g. Bohannan 1997; Clark et al. 2003), and there- 
fore could be used for selection of candidate massive stars 
for follow-up spectroscopy. Because of the huge interstellar 
extinction in the Galactic plane, the most effective channel 
for the detection of (circumstellar) shells is through imag- 
ing with modern infrared telescopes. The infrared imaging 



© 0000 RAS, MNRAS 000, 000-000 



Discovery of two new Galactic candidate luminous blue variables with WISE 3 



could also be crucial for revealing those LBV stars which 
are presently enshrouded in dusty envelopes created by ma- 
jor eruptions in the recent past and therefore are heavily 
obscured in the optical band. Application of this approach 
using the archival data of the Spitzer Space Telescope re- 
sulted in detection of dozens of new cLBV, blue supergiant, 
and Wolf-Rayet (WR) stars in the Milky Way (Gvaramadze 
et al. 2009, 2010a,b; Gvaramadze, Kniazev & Fabrika 2010c; 
Wachter et al. 2010, 2011; Mauerhan et al. 2010; Stringfel- 
low et al. 2012a,b). The majority of these shells are visible 
only in 24 /im images obtained with the Multiband Imag- 
ing Photometer for Spitzer (MIPS; Rieke et al. 2004) and 
most of them were detected in the archival data of the 24 
and 70 Micron Survey of the Inner Galactic Disc with MIPS 
(MIPSGAL; Carey et al. 2009), which mapped 278 square 
degrees of the inner Galactic plane (|6| < 1° is covered 
for 5° < / < 63° and 298° < I < 355° and |6| < 3° is 
covered for \l\ < 5°). The high angular resolution of the 
MIPS images (6 arcsec at 24 /im) even allowed us to detect 
a new circular shell in the Large Magellanic Cloud (LMC) 
and therethrough to discover a new WR star - the first- 
ever extragalactic evolved massive star detected with Spitzer 
(Gvaramadze et al. 2011a; Gvaramadze et al., in prepara- 
tion). 

Due to a limited sky coverage of the MIPSGAL and 
other Spitzer Legacy Program^], however, many evolved 
massive stars with circumstellar shells might remain unde- 
tected. The situation was somewhat improved with the re- 
cent data release of the Mid-Infrared All Sky Survey carried 
out with the Wide-field Infrared Survey Explorer (WISE; 
Wright et al. 20100. This survey provides images at four 
wavelengths: 3.4, 4.6, 12 and 22 /im, with an angular resolu- 
tion of 6.1, 6.4, 6.5 and 12.0 arcsec, respectively, and allows 
a search for infrared circumstellar shells in the regions not 
covered by Spitzer. The factor of two lower angular resolu- 
tion of the 22 /im WISE data, however, makes many of the 
shells discovered with Spitzer almost indiscernible, so that 
we expect that only the most extended yet unknown shells 
would be detected with WISE. 

In this paper, we report the discovery of two new Galac- 
tic cLBVs via detection of circular shells with WISE and 
follow-up optical spectroscopy of their central stars with 
the renewed Southern African Large Telescope (SALT). The 
new shells and their central stars are presented in Section[2] 
Section im describes our spectroscopic follow-up of the cen- 
tral stars and the data reduction. Archival and contempo- 
rary photometry of the stars is presented in Section l3~2l The 
obtained spectra are discussed in Section^ The discussion 
is given in Section[5] We summarize and conclude in Sect. [6] 



2 NEW CIRCULAR SHELLS AND THEIR 
CENTRAL STARS 

Two circular shells, which are the main subject of this pa- 
per, were discovered using the WISE 22 /im archival data. 
Both shells contain central stars (see Tabled for their co- 
ordinates), which are prominent not only in all four WISE 

1 http:/ /irsa. ipac.caltech.edu/Missions/spitzcr. html 

2 The current release of the WISE data covers 57 per cent of the 
sky. 



wavebands but also in the optical range. We call these shells 
WS1 and WS2, i.e. the WISE shells, and, in that follows, 
we will use for the central stars the same names as for the 
associated nebulae. 

FigureQ] shows the WISE, 2MASS JC s -band (Skrutskie 
et al. 2006) and DSS-II red band (McLean et al. 2000) im- 
ages of the field containing WS1. At 22 /im the shell is al- 
most circular with a radius of ~ 35 arcsec. One can also 
see an elongated (toroid-like?) region of enhanced brightness 
stretched within the shell in the northwest-southeast direc- 
tion. The 12 fim image shows the gleam of emission prob- 
ably associated with WS1 and an arc-like structure faced 
towards the centre of the 22 /im shell. The brightest part of 
this arc coincides with the region of enhanced brightness at 
the northwest of the shell. The central star of WS1 corre- 
sponds to the 2MASS J13362862-6345387 source (Skrutskie 
et al. 2006). 

WS2 was discovered during our search for bow shocks 
generated by OB stars running away from the massive 
star-forming region NGC 6357 (for motivation and the re- 
sults of this search see Gvaramadze & Bomans 2008 and 
Gvaramadze et al. 2011b, respectively). A by-product of 
the search, which used the archival data of the Midcourse 
Space Experiment (MSX) satellite (Price et al. 2001) and 
the MIPSGAL and WISE surveys, is the discovery of three 
LBV- like shells, one of which, WS2, has a central star visible 
in the optical region. Unfortunately, the MIPSGAL survey 
covers only a part of this quite extended (with a radius of 
~ 115 arcsec) shell, so that the full extent of WS2 was re- 
vealed only at 22 /im with WISE (see Fig.[2]). The MIPS 
24 /im and the WISE 22 /im images of WS2 have already 
been presented in Gvaramadze et al. (2011b). 

WS2 has a clear circular limb-brightened structure with 
a region of reduced brightness in the southeast of the shell 
and a chimney-like blow-up (more prominent at 24 /im due 
to the better angular resolution of the MIPS image) attached 
to a region of enhanced brightness in the northeast. The 
WISE 12 /im image shows the gleam of emission possibly 
related to WS2. Particularly, one can see two wings appar- 
ently extending from the central star towards the shell. One 
of them connects the star to the brightest part of the shell 
and then diverges to the "chimney", while the second (less 
prominent) wing connects the star with the west rim of the 
shell. Although it is not clear whether or not these wings 
are related to WS2, we note that similar structures can also 
be seen within the shell of the cLBV MN96 (see Fig. l(n) 
in Gvaramadze et al. 2010c or Fig. 1 in Stringfellow et al. 
2012a). 

Using the SIMBAD databas<£| and the VizieR cata- 
logue access we identified WS2 with the emission-line 
star Hen 3-1383, classified in the literature as an M3e (The 
1961), an Mep? (Sanduleak & Stephenson 1973; Allen 1978) 
and a candidate symbiotic (Henize 1976) star. Hen 3-1383 
is also known as a variable source of thermal radio emis- 
sion (Seaquist & Taylor 1990; Altenhoff, Thum & Wendker 
1994). 

The morphology of the WS1 and WS2 shells is very 
similar to that of the majority of circumstellar shells associ- 



3 http:/ /simbad. u-strasbg.fr/simbad/ 

4 http:/ /webviz.u-strasbg.fr/viz-bin/VizicR 
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Figure 1. From left to right, and from top to bottom: WISE 22, 12, 4.6 and 3.4/im, 2MASS its-band and DSS-II red band images of 
the circular shell WS1 and its central star (marked by a circle on the 2MASS and DSS-II images). The orientation and the scale of the 
images are the same. 



ated with known (c)LBVs and related evolved massive stars 
(see Gvaramadze et al. 2010c), which allowed us to expect 
that the central stars of the shells are evolved massive stars 
as well. Follow-up spectroscopy of these stars confirmed our 
expectation and leads to the discovery of two new cLBVs. 



3 OBSERVATIONS 

3.1 Spectroscopic follow-up 

Spectral observations of WS1 and WS2 were conducted with 
SALT (Buckley, Swart & Meiring 2006; O'Donoghue et al. 
2006) on 2011 June 12 and 13 during the Performance Ver- 
ification phase of the Robert Stobie Spectrograph (RSS) 
(Burgh et al. 2003; Kobulnicky et al. 2003). The long-slit 
spectroscopy mode of the RSS was used with a 1 arcsec slit 
width for observations of both stars. The grating GR900 
was used to cover the spectral range of 4200 — 7250 A with 
a final reciprocal dispersion of 0.97 A pixel -1 and FWHM 
spectral resolution of 4.74 ± 0.05 A. The seeing during the 
observations was 2.1 arcsec for WS1 and 1.6 arcsec for WS2. 
Spectra of a Xe comparison arc were obtained to calibrate 
the wavelength scale. WS1 was observed with 3 x 300 s ex- 
posures, while WS2 with 2 x 120 s and 2 x 600 s ones. Be- 
cause of significant A" s -band photometric variability of WS2 



detected in the archival databases (see Section[3]2]), we de- 
cided to re-observe this star to search for possible changes 
in its spectrum. A second spectrum of WS2 was taken on 
2011 September 10 with 3 x 800 s exposures. We used the 
same grating to cover the spectral range of 3900 - 7000 A , 
which allowed us to fill gaps in the wavelength coverage of 
the first spectrum. The seeing during this observation was 
variable, ~ 3 arcsec or worse. Spectrophotometric standard 
stars were observed after observation of each object for rel- 
ative flux calibration. 

Primary reduction of the data was done with the SALT 
science pipeline (Crawford et al. 2010). After that, the bias 
and gain corrected and mosaiced long-slit data were reduced 
in the way described in Kniazev et al. (2008). The resulting 
rich emission-line spectra of WS1 and WS2 are presented 
and discussed in Section[4] 

Equivalent widths (EWs), FWHMs and radial heliocen- 
tric velocities (RVs) of main emission lines in the spectra of 
WS1 and WS2 (measured applying the MIDAS programs; 
see Kniazev et al. 2004 for details) are summarized in Table[T] 
and Table[2] respectively. 

We note that SALT is a telescope with a variable pupil, 
so that the illuminating beam changes continuously during 
the observations. This makes the absolute flux/magnitude 
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Figure 2. From left to right, and from top to bottom: Spitzer MIPS 24 (im, WISE 22, 12, 4.6 and 3.4 /im and DSS-II red band images 
of the circular shell WS2 and its central star (marked by a circle on the DSS-II image). The MIPS image is truncated because the 
MIPSGAL survey covers only a part of the shell. The orientation and the scale of the images are the same. 



calibration impossible, even using spectrophotometric stan- 
dard stars or photometric standards. 



3.2 Photometry 

WS2 is indicated in the International Variable Star Index 
(VSX) as a suspected variable with maximum visual magni- 
tude of 12.5 (Watson 2006). Although the source of this mag- 
nitude was not indicated, we found that it comes from the 
paper by Sanduleak & Stephenson (1973), who estimated 
V — 12.5 mag as the average limiting magnitude for detec- 
tion of a continuum in their low-dispersion spectra obtained 
with an objective prism (see also Section^. Since Sanduleak 
& Stephenson (1973) did not detect the continuum in the 
spectrum of WS2, the actual visual brightness of this star 
at the moment of their observations might have been fainter 
than 12.5 mag. The variable status of WS2, however, can 
be supported by the comparison of the 2MASS and the DE- 
NIS J and K s magnitudes (see Table[3J), which shows that 
in 1998 the star brightened in the K B -band by ~ 0.84 ±0.19 
mag during ~ 2.5 months. 

To determine the contemporary photometry of WS1 
and WS2, we observed these stars with the PROMPT net- 
work, which consists of six 0.4-m robotic telescopes operat- 
ing in Chile by the University of North Carolina (Reichart 



et al. 2005). Observations were carried out on 2011 Septem- 
ber 20 and 23 with the PROMPT No. 4 telescope, which 
is equipped with BVRI filters of the Johnson-Cousins pho- 
tometric system. Transformations from the instrumental to 
the standard photometric system were obtained during sev- 
eral nights in 2010 using Landolt (1992) standard stars in 
the regions of RU 149, RU 152, and SA 98. The photomet- 
ric accuracy is ~ 0.03 mag in all the filters. The results 
are presented in Table|3j where we also give photographic 
magnitudes obtained from different sources and calibrated 
using the secondary photometric standards established from 
the PROMPT data, J and K s magnitudes from 2MASS, and 
/, J and K s magnitudes from the DENIS database (The DE- 
NIS Consortium, 2005). 

Table[3] shows that WS1 brightened in the /-band by 
0.61 ± 0.04 mag during the last 13 years. One can also see 
that this brightening was accompanied by the brightness 
increase in the i?-band by 0.68 ± 0.10 mag. Thus, the star 
apparently became redder during the last 13 — 18 years. 
Similarly, Table[3] shows that WS2 became fainter in the V- 
band by 0.47 ±0.05 mag during 10 days of 2011 September, 
and brightened in the /-band by 0.41 ± 0.05 mag during 
three months since the first spectrum was taken, i.e. the 
star appears to have significantly reddened on a time-scale 
of several months. Table[3]also shows that WS2 experienced 
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Table 2. Equivalent widths (EWs), FWHMs and radial heliocentric velocities (RVs) of some detected lines in the spectra of WS2 taken 
on 2011 June 13 (1) and 2011 September 10 (2). RVs of lines noticeably affected by P Cygni absorptions are starred. 



Arj(A) Ion 


EW(A) 


FWHM(A) 


RV 




(A) 




(A) 


(krns- 1 ) 




(i) 


(2) 










4340 H7 


6.92±0.35 


8.88±0.49 


6.76±0.40 


6.23±0.40 


129±11* 


107±12* 


4471 He I 


3.63±0.21 


2.00±0.07 


7.01±0.47 


6.61±0.31 


134±13* 


76±8* 


4554 Si III 


2.66±0.10 


1.52±0.28 


9.23±0.27 


8.13±1.96 


54±8 


44±50 


4631 N n 


1.41±0.05 


0.66±0.15 


8.08±0.31 


7.97±0.49 


76±7 


63±10 


4713 He I 


1.98±0.05 


2.20±0.09 


7.22±0.13 


7.69±0.35 


96±3 


85±9 


4861 H/3 32.69±0.56 


31.61±0.49 


7.08±0.14 


7.03±0.12 


72±4* 


58±3* 


4922 He I 


4.22±0.38 




6.21±0.64 


- 


164±17* 


- 


5016 He I 


7.90±0.45 


10.71±0.58 


8.34±0.56 


8.32±0.52 


176±14* 


179±13* 


5128 Fe in 


1.36±0.10 


1.45±0.11 


5.72±0.42 


6.15±0.47 


55±11* 


77±12* 


5156 Fc in 


1.72±0.07 


2.12±0.10 


6.11±0.23 


6.52±0.33 


75±6* 


91±8* 


5171 Fc n 


3.25±0.09 


3.55±0.10 


7.64±0.21 


8.29±0.26 


-36±5* 


-26±6* 


5319 Fc n 


2.41±0.10 


2.61±0.09 


8.55±0.20 


9.82±0.54 


-56±5 


-65±8 


5334 [Fc n] 


0.70±0.11 


0.65±0.09 


11.63±0.90 


9.89±1.44 


— 24±18 


-40±32 


5376 [Fe n] 


0.58±0.09 


0.64±0.06 


13.56±0.32 


11.25±0.64 


-16±7 


7±15 


5755 [N n] 


2.07±0.09 


1.98±0.08 


13.11±0.40 


11.24±0.48 


-122±18 


-117±10 


5834 Fe III 


0.72±0.08 


0.59±0.19 


8.72±0.23 


4.81±1.69 


21±5* 


-11±37* 


5876 He I 13.04±0.26 


16.42±0.35 


6.75±0.15 


7.06±0.17 


30±3* 


42±4* 


5979 Fe m 


1.70±0.09 




10.93±0.23 




-53±5 




6000 Fe in 


0.29±0.06 




7.51±0.65 




-92±13 




6033 Fe m 


0.54±0.05 




6.96±0.15 




-65±3 




6347 Si n 


1.43±0.06 


1.29±0.05 


9.73±0.27 


10.37±0.30 


13±5* 


42±5* 


6371 Si n 


0.68±0.05 


0.66±0.04 


7.45±0.12 


7.49±0.22 


-38±2 


-1±4 


6482 N II 


1.48±0.20 


1.50±0.20 


9.26±1.34 


9.49±1.39 


188±26* 


201±27* 


6563 Ha 136.00±2.55 108.98±2.10 


11.52±0.25 


13.73±0.30 


43±5* 


55±6* 


6678 He I 


7.83±0.45 


8.37±0.46 


7.77±0.51 


7.96±0.50 


54±10* 


68±10* 


7065 He I 12.28±0.44 




9.17±0.38 




11±7* 




7155 [Fe n] 


1.92±0.16 




16.82±0.86 




-64±15 




Archival and contemporary photometry of WS1 and WS2. 


WS1 


Date 


B 


V 


R 


I 


J 




1976 March 8^ 


17.50±0.10 












1979 June 7^ 








12.60±0.10 






1984 February 22* 1 ) 






14.10±0.10 








1993 June 20^ 






14.30±0.10 








1998 July 1< 2 ) 








12.79±0.03 


10.56±0.10 


8.90±0.11 


2000 February 22< 3 ) 










10.33±0.02 


8.90±0.02 


2011 June 12( 4 > 




15.25±0.04 










2011 September 23^ 


17.36±0.03 


15.31±0.03 


13.62±0.03 


12.18±0.03 






WS2 


Date 


B 


V 


R 


I 


J 


K s 


1958 April 17' 1 ) 


17.70±0.10 












1976 May 3lW 


16.90±0.10 












1980 May lM 








10.20±0.10 






1984 April 9^) 






12.20±0.10 








1991 August 5W 






12.70±0.10 








1998 August 10( 3 > 










6.62±0.04 


4.70±0.02 


1998 October 29< 2 > 








9.86±0.05 


6.18±0.12 


3.86±0.19 


2011 June 13( 4 > 








10.13±0.04 






2011 September 10< 4 ) 




14.10±0.04 










2011 September 20< 5 > 


17.80±0.03 


14.57±0.03 


11.89±0.03 


9.72±0.03 






(1) USNO B-l (Monet et al. 2003); 


(2) DENIS; (3) 2MASS; (4) SALT; (5) PROMPT. 
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Table 1. Equivalent widths (EWs), FWHMs and radial helio- 
centric velocities (RVs) of some detected lines in the spectrum of 
WS1 (taken on 2011 June 12). RVs of lines noticeably affected by 
P Cygni absorptions are starred. 



A (A) Ion EW(A) FWHM(A) RV 

(A) (A) (kms" 1 ) 



4340 


H7 


4.7±0.2 


4.97±0.21 


83±6* 


4471 


He I 


3.4±0.3 


8.08±0.79 


205±22* 


4554 


Si in 


0.5±0.1 


5.18±0.72 


198±20* 


4631 


N n 


0.2±0.1 


6.16±0.22 


102±5 


4658 


[Fc in] 


0.12±0.04 


6.29±0.42 


-103±8 


4686 


He II 


0.4±0.1 


10.82±0.43 


94±8 


4713 


He I 


1.1±0.1 


5.82±0.19 


90±5 


4861 


H/3 


16.3±0.3 


6.53±0.12 


42±3* 


4922 


tic I 


1.3i0.2 


a KO_l_n o 
4.0Z±0.DO 




5016 


He I 


4.3±0.2 


5.98±0.27 


79±7* 


5128 


Fc in 


0.4±0.1 


4.25±1.20 


139±30* 


5156 


Fc in 


0.2±0.1 


3.33±1.59 


136±39* 


5755 


[N ii] 


0.8±0.2 


8.14±0.72 


-41±16 


5834 


Fe ill 


0.5±0.2 


6.65±0.37 


-1±8 


5876 


He I 


16.5±0.5 


6.29±0.19 


64±4* 


5979 


Fc in 


0.7±0.2 


6.63±0.21 


-16±4 


6000 


Fc in 


0.2±0.1 


4.06±0.54 


-6±12 


6033 


Fc in 


0.7±0.2 


6.27±0.10 


-14±2 


6347 


Si ii 


1.3±0.2 


12.37±0.49 


-56±9 


6371 


Si ii 


0.5±0.1 


6.84±0.23 


-31±5 


6482 


N ii 


1.0±0.1 


5.98±0.24 


76±5* 


6563 


Ho 


63.2±0.9 


7.86±0.12 


32±2* 


6611 


N ii 


0.7±0.1 


9.13±0.49 


-42±8 


6678 


He I 


6.9±0.3 


6.27±0.28 


69±5* 


7065 


He I 


10.7±0.3 


7.20±0.17 


30±3* 



Table 4. Details of central stars associated with two circular 
shells discovered with WISE. 



WS1 WS2 (= Hen 3-1383) 



RA(J2000) 


13 h 36 m 28?63 


17 h 20 m 31?76 


Dec.(J2000) 


-63°45'38'/7 


-33°09'49'/0 


I 


307.8856 


353.5275 


b 


-1.3230 


2.2005 


B (mag) 


17.36±0.03 


17.80±0.03 


V (mag) 


15.31±0.03 


14.57±0.03 


[3.4] (mag) 


8.03±0.03 


4.18±0.05 


[4.6] (mag) 


7.59±0.02 


2.79±0.03 


[12] (mag) 


6.76±0.03 


2.58±0.02 


[22] (mag) 


4.62±0.04 


1.65±0.03 



"PROMPT; 6 2MASS; C WISE. 



significant (0.5 — 0.9 mag) variability in the B- and 7?-bands 
on time-scales of decades, although the poor time cadence 
leaves the possibility that the brightness changes in these 
bands occurred on much shorter time-scales. 

The details of WS1 and WS2 are summarized in Tabled 
where we give their 2MASS coordinates, the contemporary 
B and V magnitudes, and, for the sake of completeness, 
the WISE magnitudes from the WISE Preliminary Release 
Source Catalogj. 



4 WS1 AND WS2 - NEW GALACTIC 
(CANDIDATE) LBV STARS 

Figures[3E] present the normalised spectra of WS1 and WS2 
(taken on 2011 June 12 and 13, respectively) in the blue, 
green and red regions, where most of the identified lines 
are marked. For the sake of comparison, we also show in 
these figures the corresponding spectra of the prototype 
LBV P Cygni (taken from Stahl et al. 1993 and degraded 
to the resolution of our spectrajf] and the recently discov- 
ered Galactic cLBV MN112 (Gvaramadze et al. 2010b). The 
spectra are arranged in such a way that objects with the sup- 
posedly cooler effective temperature (see below) are placed 
closer to the top of the figures. One can see that all four spec- 
tra are very similar, which suggests the LBV classification 
for WS1 and WS2 as well. To identify the lines indicated in 
the figures, we mostly used the spectral atlas of P Cygni by 
Stahl et al. (1993). 

All four spectra are dominated by strong emission lines 
of H and He I; some of them show P Cygni profiles. Both Ha 
and H/3 have prominent wings. Other emission lines present 
in all four spectra are numerous metal lines of N II, Fein and 
Sin, and the carbon doublets at AA 6578, 6583 and AA 7231, 
7236. The dichotomy of Fein lines observed in the spectra 
of P Cygni (Stahl et al. 1993) and MN112 (Gvaramadze et 
al. 2010b) is present in the spectra of WS1 and WS2 as well: 
lines with low multiplet numbers (e.g. Fern AA5127, 5156) 
have distinct P Cygni profiles, while those with higher mul- 
tiplet numbers (e.g. Fe ill AA5920 — 6033) are purely in emis- 
sion. Some of the N II lines (e.g. AA6151, 6168) are also purely 
in emission. The only forbidden line common for all four 
spectra is that of [Nil] A5755, although the presence of weak 
[Nil] AA6548, 6583 lines in the Ha emission wings cannot be 
excluded. 

Numerous diffuse interstellar bands (DIBs) are present 
in the spectra. In the blue region they are AA4428, 4726 and 
4762, in the green and red the strongest DIBs are at 5780, 
5797, 6379 and 6613 A. Strong absorptions visible at A > 
6800 A are telluric. The sodium doublet in all four spectra is 
of interstellar origin. In addition to the interstellar lines, the 
spectrum of WS2 (as well as that of P Cygni; Stahl et al. 
1993) shows blueshifted (Aw ~ —220 kms -1 ) components 
of circumstellar origin. 

There are, however, several important distinctions in 
the spectra. The blue spectrum of WS1 shows the presence 
of a weak emission of Hell A4686 and the high-ionization 
forbidden line [Fein] A4658, which indicates that this star 
is somewhat hotter than the other three stars. The presence 
of these lines makes the spectrum of WS1 almost identical 
to that of the bona fide LBV AG Car during the epoch of a 
minimum in 1985-1990, when it was a WN11 star (Stahl et 
al. 2001) with an effective temperature of ~ 22000 - 23000 
K (Groh et al. 2009a). Thus, WS1 can be classified as a 
WN11 star. The WN11 classification of WS1 also follows 
from the position of this star on the diagrams showing a 
comparison of EWs of He I A5876 versus Hell A4686 lines 
and Hen A4686 EW versus FWHM for late WN-type (WNL) 
stars in the Milky Way and the LMC (see Fig. 1 of Crowther 
& Smith 1997). With EW(5876)=16.5 A, EW(4686)=0.4 A, 
and FWHM(4686) = 10.8 A (see TableHJ, WS1 lies in the 



5 http:/ /irsa. ipac.caltech.edu/cgi-bin/Gator/nph-dd 6 http:/ /www. lsw.uni-hcidclbcrg.de/users/ostahl/pcyg/ 
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Figure 3. Comparison of normalised blue spectra of WS1 and WS2 with those of the prototype LBV P Cygni and cLBV MN112. The 
principal lines and most prominent DIBs are identified. 



region occupied by WN11 stars. We conclude, therefore, that 
WS1 is a cLBV near a brightness minimum. 

The spectrum of WS2 shows several permitted and for- 
bidden emission lines of singly-ionized iron, which are ab- 
sent or too weak to be detected in the spectra of other 
three stars. Of these lines the most prominent are Fell 
AA5171, 5275, 5319, 5899, 6382 and [Fen] AA5334, 5376, 7155. 
Some of the Fell lines (AA5892 - 5897) are blended with the 
blueshifted component of the Nal Dl line, while some oth- 
ers (e.g. AA4922, 5016) are blended with the He I lines. The 
latter results in an anomalously high intensity ratio of He I 
A4922 to He I A4471 compared to the other three stars (cf. 
Hillier et al. 1998). Several emissions of Alii are also de- 
tected in the spectrum, of which the most prominent are at 
AA6226 - 6243 and A7042. Overall, the spectrum of WS2 is 
very similar to those of AG Car in 2003 January (see Fig. 3 
in Groh et al. 2009a) and the cLBV HD 316285 (Hillier et 
al. 1998). Since the singly-ionized iron lines appeared in the 
spectrum of AG Car two years after a visual minimum in 
2001 (when the effective temperature of AG Car decreased 
from ~ 17000 to 14000 K; Groh et al. 2009a), one can con- 
clude that WS2 is also near the brightness minimum and 
currently is either at the onset of the rise to the maximum 
or on the way to the minimum. 

Figures r3l5l show that the emission lines of WS2 are 
much broader than those of WS1 (cf. also Tables [1] and [2} 



and the two other stars, which might be caused by a higher 
wind velocity, Voo, of WS2. Indeed, this inference is sup- 
ported by measurements of the FWHM of the [N n] A5755 
line, which can be used as a measure of Voa (e.g. Crowther, 
Hillier & Smith 1995). After correction for instrumental 
width (~ 4.74 ± 0.05 A), we found ~ 640 ± 20 kms -1 
for WS2 and ~ 350 ± 40 km s" 1 for WS1. A similar value of 
Voo for WS2 also follows from the width of the flat-topped 
lines [Fell] AA 5334, 5376, and 7155: the mean value of Doo 
derived from these three lines is ~ 690 ± 15 kms -1 . For 
comparison, the wind velocities of P Cygni and MN112 are 
~ 200 kms -1 (e.g. Barlow et al. 1994; Najarro, Hillier & 
Stahl 1997) and ~ 400 kms" 1 (Gvaramadze et al. 2010b), 
respectively. Although the high wind velocity inferred for 
WS2 is more typical of WN8-9h stars (Hamann, Grafener 
& Liermann 2006), we note that similar velocities were also 
measured for the Homunculus nebula of r\ Car (e.g. Allen 
& Hillier 1993; Meaburn, Walsh & Wolstencroft 1993) and 
for the LBV-like component of the WR binary HD 5980 in 
the Small Magellanic Cloud during its slow wind eruptive 
phase in 1994, when this component was transformed from 
a WR star in a LBV (Koenigsberger et al. 1998). This trans- 
formation was accompanied by drastic changes of the wind 
velocity, which reduced from ~ 3000 to 500 kms -1 and af- 
ter the eruption increased to the initial value. Thus, it is 
possible that WS2 has either experienced a giant eruption 
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Figure 4. Comparison of normalised green spectra of WS1 and WS2 along with spectra of the prototype LBV PCygni and cLBV 
MN112. The principal lines and most prominent DIBs are identified. 



in the past, so that the inferred high velocity corresponds 
to the polar outflow a la the Homunculus nebula, or that it 
is caught shortly before or soon after an LBV-like eruption, 
similar to that occurred in HD 5980 in 1994. In the latter 
case, one can expect that WS2 will show major changes in 
its spectrum in the near future. 

To search for possible spectral variability of WS2, we 
obtained a second spectrum of this star ~ 3 months after 
the first observation. Table[2] shows that although most lines 
remain invariable within their margins of error, some of them 
have changed their characteristics significantly. The most 
noticeable changes occurred in the EWs of the He I AA 5016, 
5876 lines (which became larger by ~ 30 — 40 per cent) and 
in the EW and the FWHM of the Ha line (which decreased 
and increased by ~ 20 per cent, respectively). To illustrate 
these changes graphically, we divided the first normalised 
spectrum of WS2 by the second one. Figs.|6]and[7]show the 
result for two spectral ranges, 5700 — 6000 A and 6300 — 
6750 A, respectively. Changes in EWs, RVs and the line 
profiles manifest themselves in the origin of artefacts above 
and below continuum. Some of these artefacts show P Cygni- 
like profiles (e.g. at the positions of the Ha and He I A6678 
lines), while others (e.g. He I A5876, 6678) have an additional 
inverse P Cygni-like profile. 

Tabled also shows that the FWHMs of the [Nil] A5755 
and the flat-toped [Fell] lines are systematically reduced, 



which indicates the reduction of the wind velocity as well. 
Using the FWHMs of these lines, we derived the mean wind 
velocity of ~ 540 ± 20 kms -1 , which implies that Vao de- 
creased by ~ 100 ± 20 kms" 1 during the last three months 
(recall that this velocity decrease is accompanied by the sig- 
nificant reddening of the star). 

To conclude, the similarity between the spectra of WS1 
and WS2 and those of the well-known LBVs P Cygni and 
AG Car, along with the significant short-term photomet- 
ric variability and the presence of the circumstellar shells 
strongly argue for the LBV classification of both stars. For 
WS2 this classification is further supported by the changes 
in the spectrum. We hope that spectroscopic and photomet- 
ric monitoring of WS1 and WS2 will allow us to confirm the 
LBV status of these stars in the foreseeable future. 



5 DISCUSSION 

5.1 Classification of WS2 

We reported the detection of two circular shells with WISE 
and the results of follow-up optical spectroscopic observa- 
tions of their central stars with SALT, which lead to the 
discovery of two new Galactic cLBVs, WS1 and WS2. WS2 
was known long ago as an emission-line star Hen 3-1383 and 
was classified in the literature as an M-type star, which is 
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Figure 5. Comparison of normalised red spectra of WS1 and WS2 along with spectra of the prototype LBV P Cygni and cLBV MN112. 
The principal lines and most prominent DIBs arc identified. 
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Figure 6. The ratio of the normalised (green) spectra of WS2 
obtained within ~ 3 months from each other. The positions of 
principal lines and most prominent DIBs are indicated. 



Figure 7. The ratio of the normalised (red) spectra of WS2 ob- 
tained within ~ 3 months from each other. The positions of prin- 
cipal lines and most prominent DIBs arc indicated. 



completely discrepant with our classification of this star as a 
cLBV. Let us discuss the possible origin of the discrepancy. 

The existing spectral classifications of WS2 were based 
on the observations carried out ~ 30 — 60 years ago, so that 



one can wonder whether WS2 has evolved blueward dur- 
ing the last several decades? Although there are no known 
examples of such a radical change in the spectral type of 
any know star, one cannot completely exclude this possibil- 
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ity. For instance, in the recent past WS2 might be a yel- 
low hypergiant star, which, like the M33's hypergiant star 
Var A (Humphreys, Jones & Gehrz 1987), experienced peri- 
odic changes of spectral type from pseudo M-type to G-type 
during excursions on the cool side of the Yellow Void (de 
Jager & Nieuwenhuijzen 1997). These excursions were ac- 
companied by episodes of enhanced mass loss due to atmo- 
spheric instability (de Jager 1998). Finally, several decades 
ago WS2 made a last blueward swing and after crossing the 
Yellow Void found themselves close to the S Doradus insta- 
bility strip (Wolf 1989), where we currently observe it as an 
LBV. WS2 therefore might be a more evolved analog of the 
blueward evolving yellow hypergiant IRC+10 420 (Oudmai- 
jer 1998), which most probably will turn into an LBV in 
the close future (Humphreys, Davidson & Smith 2002). The 
characteristic timescale of formation and dispersal of obscur- 
ing envelopes around yellow hypergiants ranges from years 
to decades, which is comparable to the time period elapsed 
between our and the previous spectroscopic observations of 
WS2. 

The above considerations imply that WS2 has emerged 
as an LBV only a few decades ago. This implication, how- 
ever, is inconsistent with the presence of the extended 
(parsec-scale; see below) circumstellar shell around this star 
because shells of this scale are rather typical of LBV and 
WNL stars, but were never observed around known yellow 
hypergiants (e.g. Humphreys et al. 1997; Jura & Werner 
1999; Smith et al. 2001; Lagadec et al. 2011). Although one 
cannot exclude the possibility that WS2 became an LBV 
star relatively long ago and that in the recent past it ex- 
perienced a giant eruption leading to the formation of an 
extended pseudo-photosphere, which pushed the star on the 
cool side of the Yellow Void (cf. Smith et al. 2004) and which 
dispersed during the last several decades, we believe that the 
origin of the discrepancy between the former and the present 
spectral classifications of WS2 is rather more prosaic. Let us 
discuss the sources of the previous classifications. 

WS2 was identified as an emission-line star by The 
(1961), who found a strong Ha emission in its spectrum 
(obtained in 1960). The (1961) classified the star as M3 and 
suggested that "it is possibly a long-period variable star", 
although no information on the basis of these findings was 
provided. The visual brightness of WS2 was measured to be 
14 mag. 

Sanduleak & Stephenson (1973) classified WS2 as a 
Mep? star on the basis of detection of an extremely strong 
Ha emission (with no evidence of H/3 and He n A4686 emis- 
sions) and the suspected presence of TiO bands. WS2 does 
not show evidence of a continuum on the survey plates 
(taken in 1967-1968), from which Sanduleak & Stephenson 
(1973) concluded that the star was fainter than 12.5 mag at 
that time. This magnitude is quoted in several subsequent 
catalogues (e.g. in Watson 2006, where the star is indicated 
as a possible symbiotic variable of the Z Andromedae type). 

The TiO bands (and the continuum as well) were not 
visible in the spectrum taken in 1949-1951 by Henize (1976), 
who interpreted the presence of the very strong Ha emission 
as the indication that WS2 might be a symbiotic star. 

The spectral classifications of WS2 given in the all 
three above-mentioned papers were based on objective- 
prism spectra. The first (low-dispersion) slit spectrum of 
WS2 was obtained in 1977 by Allen (1978), who suspected 



the presence of weak TiO absorption and tentatively clas- 
sified the star as a peculiar M-type emission-line star (i.e. 
Mep?), because the Hi and He I emissions in its spectrum 
"are very strong for an Me star". 

To summarize, the existing historical spectroscopic data 
show that WS2 is an emission-line star with the very strong 
Ha emission (detected in all spectra starting from 1949 till 
1977) and the Hi emission (detected in 1977 in the first slit 
spectrum of the star). TiO bands were suspected in 1967- 
1968 and 1978, and were not visible in 1949-1951. Although 
with this information we cannot exclude the possibility that 
WS2 was an pseudo M-type emission-line star ~ 30 — 60 
years ago, it is plausible that the previous spectral clas- 
sifications are erroneous because of the low quality of the 
spectroscopic material, aggravated by high visual extinction 
(Av ~ 10 mag; see below) towards the star. 

5.2 Reddening towards WS1 and WS2 

We now turn to estimate the reddening towards WS1 and 
WS2. First, we matched the dereddened spectral slopes of 
these stars with those of stars of similar effective temper- 
ature. Using the Stellar Spectral Flux Library by Pickles 
(1998) and assuming the effective temperatures of WS1 
and WS2 of 22000 K and 14000 K (see Section©, respec- 
tively, we found the colour excess E(B — V) = 2.40 mag 
for WS1 and E(B - V) = 3.37 mag for WS2. These esti- 
mates only slightly depend on the assumed effective tem- 
perature. By varying the temperature in a quite wide range 
(~ 12000 - 26000 K for WS1 and ~ 11000 - 23000 K for 
WS2), we found E(B - V) ~ 2.37 ± 0.10 mag for WS1 and 
E(B-V) ~ 3.36 ±0.10 mag for WS2. We also compared the 
spectral slopes of WS1 and WS2 with that of AG Car (using 
the spectra of AG Car obtained by us on 2009 June 1 and 
2011 May 19), and found that they nicely match each other 
if one dereddens the spectra of WS1 and WS2, respectively, 
by 1.53 - 1.66 and 2.61 - 2.73 mag in E(B - V). Adopt- 
ing E(B — V) = 0.63 mag for AG Car (Humphreys et al. 
1989), we found E(B -V) = 2.16 - 2.29 mag for WS1 and 
E(B -V) = 3.24 - 3.36 mag for WS2, which are fairly well 
agree with the above estimates. 

The reddening can also be estimated by using the corre- 
lation between the intensity of the DIBs and E(B — V) (see 
Herbig 1995 for a review). In Table[5]we provide estimates of 
E(B-V) for both stars based on EWs of DIBs at A5780 and 
A5797, and the relationships given in Herbig (1993). We did 
not find this approach useful because of the large spread in 
the resulting estimates of E(B — V) and their inconsistency 
with those derived from the spectral slopes. The discrep- 
ancy could be caused by foreground regions of enhanced 
number density of carriers of the DIBs. Another possibil- 
ity is that the circumstellar material around WS1 and WS2 
produces the DIB carriers and therethrough contributes to 
the strength of the DIBfl (cf. Tug & Schmidt-Kaler 1981; 
Le Bertre & Lequeux 1993; Heydari-Malayeri et al. 1993). 
In what follows we assume E(B — V) = 2.4 mag for WS1 
and E(B — V) = 3.4 mag for WS2, and use these figures 

7 We measured the EWs of DIBs in both spectra of WS2 and 
found that they agree with each other within the margin of errors 
(see also Fig. [6}. 
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Table 5. Estimates of the colour excess, E(B — V), towards WS1 
and WS2 based on the equivalent widths (EWs) of diffuse inter- 
stellar bands (DIBs) at A5780 and A5797 in the spectra of these 
stars. 





WS1 




WS2 




DIB 

(A) 


EW(A) 
(A) 


E(B - V) 
(mag) 


EW(A) 
(A) 


E(B-V) 
(mag) 


5780 
5797 


1.79±0.22 
0.83±0.19 


3.38±0.52 
6.06±0.88 


1.78±0.12 
0.91±0.12 


3.36±0.46 
6.69±0.85 



to constrain the distances to and the luminosities of these 
stars. Using the standard ratio of total to selective extinc- 
tion Rv = 3.1, we obtained Ay — 7.4 mag for WS1 and 
A v ~ 10.5 mag for WS2. 

5.3 Location of WS1 and WS2 

To estimate the distance to WS1 we use the empirical fact 
that LBVs in the hot state are located on the S Doradus 
instability strip (Wolf 1989). Assuming that the current ef- 
fective temperature of WS1 is ~ 22000 - 23000 K (see Sec- 
tion[4]), we found the minimum luminosity of this star of 
\og(L / Lq) ~ 5.8. Then assuming that the bolometric cor- 
rection of WS1 is equal to that of AG Car during the epoch 
of minimum in 1985-1990, -2.4 mag (Groh et al. 2009a), we 
found the absolute visual magnitude My ~ —7.3 mag and 
the distance of ~ 11 kpc (which places WS1 on the far side 
of the Carina-Sagittarius arm). At this distance the linear 
radius of the shell is ~ 1.8 pc (i.e. a figure typical of LBV 
shells; e.g. Weis 2001), while the distance from the Galactic 
plane is ~ 240 pc, which implies that WS1 is a runaway star 
(e.g. Blaauw 1993). The luminosity of WS1 would be higher 
if this star left the S Doradus instability strip and now un- 
dergoes an outburst or recovers after it. Correspondingly, 
the distance to the star would be larger, while its runaway 
status would be even more unavoidable. 

In Section^ we mentioned that WS2 was detected not 
far (~ 1.4 degree to the northwest) from the star-forming re- 
gion NGC 6357 and its associated massive star clusters Pis- 
mis 24 and AH03 J1725— 34.4, which are located on the close 
side of the Carina-Sagittarius arm at a distance of ~ 1.7 kpc 
(Gvaramadze et al. 2011b). At the same time, WS2 is pro- 
jected on the eastern periphery of the Sco OB4 association 
(see Fig. [8]), which is located at a distance of ~ 1.1 kpc 
(Kharchenko et al. 2005). Let us discuss whether or not one 
of these stellar systems might be the birthplace of WS2. 

The possible association of the WC7 star WR93 with 
Pismis 24 (Massey et al. 2001) implies that the age of this 
cluster is > 2 — 2.5 Myr, so that the cluster is old enough 
for some of its (most massive) members to enter the LBV 
phase. On the other hand, the ~ 40 pc separation (in projec- 
tion) of WS2 from Pismis 24 would imply that WS2 is a run- 
away star. Assuming that WS2 was ejected from Pismis 24 
and that the ejection event occurred ~ 2 Myr ago, one has 
the peculiar transverse velocity of this star of ~ 20 kms -1 . 
The peculiar velocity of WS2, however, would be several 
times higher if the star was ejected less than 2 Myr ago 




Figure 8. MSX 8.3 fim image of the star-forming re- 
gion NGC 6357 (containing the star clusters Pismis 24 and 
AH03 J1725— 34.4; indicated by diamonds) and its environments. 
The position of WS2 is marked by a small circle. The approxi- 
mate boundary of the Sco OB4 association is shown by a large 
dashed circle. The Galactic plane is shown by a dashed line. 

and if its space velocity has a significant radial component. 
Such peculiar velocities can easily be achieved in the course 
of three-body dynamical encounters between very massive 
binary and single stars in the dense cores of young star 
clusters (e.g. Gvaramadze & Gualandris 2011). Detection 
of numerous very massive runaway stars around very young 
(~ 1 — 2 Myr) massive star clusters (Evans et al. 2010; Gvara- 
madze, Kroupa & Pflamm-Altenburg 2010d; Bestenlehner et 
al. 2011; Roman-Lopes, Barba & Morrell 2011) supports the 
theoretical expectation that the dynamical evolution of such 
clusters is dominated by the most massive stars, so that the 
ejected stars are also preferentially the most massive ones. 
Unfortunately, the existing proper motion measurements for 
WS2 are very unreliable, so that we did not attempt to cal- 
culate its peculiar transverse velocity and thus to check its 
runaway status. 

Assuming that the bolometric correction of WS2 is 
equal to that of AG Car during the epoch of minimum in 
2003, -1.2 mag (Groh et al. 2009a), and adopting V = 14.1 
mag and the distance of 1.7 kpc, we found My ~ —7.6 
mag and log(L/ Lq) ~ 5.43- These estimates imply that 
WS2 might belong to a group of low-luminosity LBVs with 
log(Z//Z/Q) ~ 5.2 — 5.6 (Humphreys & Davidson 1994), which 
in turn would imply that the progenitor of WS2 has evolved 
through a red supergiant stage (Humphreys & Davidson 
1994; Smith et al. 2004). The latter implication, however, 
means that the initial mass of WS2 was < 25 — 30 Mq , and 
that its age is > 5 Myr, which is inconsistent with the young 
age of Pismis 24, and argues against the association between 
these two objects. 

If WS2 is a member of the Sco OB4 association (d = 1.1 



8 Note that these estimates should be considered approximate 
because of the significant variability of WS2 (see Tablcl3t . 
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kpc; Kharchenko et al. 2005), then Mv — —6.6 mag and 
log(L/L©) ~ 5.1. The resulting luminosity is therefore well 
below the minimum luminosity found for any known (c)LBV 
and would place WS2 on the hot side of the S Doradus in- 
stability strip. Therefore, WS2 is most likely a background 
star projected on Sco OB4 by chance. 

It is also possible that WS2 is located in one of the next 
two arms out, the Crux-Scutum Arm or the Norma Arm. In 
this case, the distance to WS2 is either ~ 3 or 5 kpc, which 
correspond to Mv ^ —8.8 mag and \og(L/L & ) ~ 5.9 or 
Mv — —9.9 mag and \ag(L/L®) ~ 6.4, respectively. Note 
that for P Cygni and AG Car (during the minimum phase) 
these parameters are equal, respectively, to —8.0 mag and 
5.8 (Najarro et al. 1997) and -8.1 mag and 6.2 (Groh et al. 
2009a). At the distance of 3 — 5 kpc the linear radius of the 
shell is ~ 1.6 — 2.7 pc, which is still within the range of radii 
derived for LBV shells. Placing WS2 at a larger distance, 
say at 8 kpc, increases its Mv and \og(L / Lq) to ~ —10.9 
mag and ~ 6.8, respectively. This would make WS2 the most 
luminous (c)LBV known to date, while its separation from 
the Galactic plane and the radius of the shell would be ~ 
300 pc and ~ 4.4 pc, respectively. Taken together, these 
extreme figures make the large distances to WS2 unlikely. 
To conclude, we deem that the distance to WS2 is either 
~ 3 or 5 kpc, although we cannot exclude the possibility 
that WS2 was expelled from a star cluster (of appropriate 
age) hidden in the molecular cloud associated with the star- 
forming region NGC 6357 (cf. Gvaramadze et al. 2011b). 

5.4 WS1 and WS2 as runaway stars 

We now discuss the possibility that the LBV activity of WS1 
and WS2 and their possible runaway status might be related. 

The fraction of runaway O stars is found observationally 
to be about 10-30 per cent (Gies 1987; Stone 1991; Zinnecker 
& Yorke 2007). Gies (2007) summarizes the evidence for 
runaway O stars to often be rapid rotators. This idea is 
enhanced by the finding of Walborn et al. (2010, 2011) that 
most of the spectroscopically identified fast-rotating O stars, 
the On stars, are runaways or are located in the field (and 
therefore most probably are runaways as well; Gvaramadze 
et al. 2010d; Weidner et al. 2011). 

Theoretically, a connection between the runaway sta- 
tus and fast rotation of massive stars is expected. Mass 
transfer in a close binary may spin up the mass receiver, 
which becomes a runaway star upon the supernova explo- 
sion of the mass donor (Petrovic, Langer & van der Hucht 
2005; Eldridge, Langer & Tout 2011). Alternatively, rapid 
rotators could also be produced in the course of dynamical 
few-body encounters in young massive star clusters, either 
because of merging of two (or several) stars or due to the 
strong tidal interaction (Alexander & Kumar 2001). In both 
cases, the fast-rotating star could be ejected from the cluster 
(e.g. Leonard 1995; Vanbeveren et al. 2009). Moreover, dy- 
namical encounters can also produce runaway binaries (e.g. 
Leonard & Duncan 1990; Kroupa 1998), which then can pro- 
duce rapid rotators by mass transfer or merging. 

While the physical mechanism for LBV outbursts is still 
unclear, there is consensus that the Eddington limit plays an 
important role in these events (Lamers & Fitzpatrick 1988; 
Smith & Conti 2008) and that the critical rotational veloc- 
ity of stars near their Eddington limit becomes very small 



(Langer 1997, 1998). The latter fact has two important con- 
sequences. First, it implies that rotation must be important 
in LBV outbursts (Langer 1997; Langer et al. 1999; Maeder 
& Meynet 2000), which is supported by the observation (e.g. 
Nota et al. 1995; Weis 2011) that the youngest LBV nebu- 
lae have a bipolar shape (note that such nebulae may be- 
come more spherical with time). Second, while fast rotation 
may not be a sufficient condition to initiate the LBV phe- 
nomenon, rapid rotators are more likely to become LBVfl 
In stars of relatively low initial mass, fast rotation may even 
trigger the LBV activity, i.e. such stars would not become 
LBVs if they were slow rotators. 

In this connection, we note that tidal stripping of en- 
velops of evolved massive stars in close dynamical encounters 
might drive stars towards the Eddington limit and thereby 
to trigger their LBV activity. This would allow stars at the 
low- mass end of the range of massive stars to manifest them- 
selves as LBVs, even if they cannot approach the Eddington 
limit in the course of an unperturbed evolution (cf. Smith 
et al. 2011). Some of these stars might also be ejected in the 
field. 

Recently, Groh, Hillier & Damineli (2006) and Groh 
et al. (2009b) found that the bona fide LBVs AG Car and 
HRCar appear to rotate close to their break-up velocities. 
Groh et al. (2009b) suggested that the fast rotation is typi- 
cal for all LBVs with strong variability on short (~ decades) 
timescales, while the dormant or ex-LBVs (e.g. P Cygni 
and HD 168625) rotate more slowly. These slowly-rotating 
(dormant /ex-) LBVs might have been fast rotators in the 
past, which have lost a significant fraction of their angular 
momentum due to one or several episodes of violent mass 
ejection. 

When rapid rotators are more likely to become LBVs, 
and rapid rotators are often runaway stars, we conclude that 
LBVs may often be runaway stars. This is supported by 
the observation that AG Car and HRCar, like many other 
(c)LBVs, are located in the field, i.e. far from known star 
clusters, and therefore are likely runaways. We thus suggest 
that the LBV activity of WS1 and WS2 might be directly 
related to their runaway status. 



6 SUMMARY AND CONCLUSION 

We have identified two new Galactic candidate LBVs via de- 
tection of mid-infrared circular shells with WISE and follow- 
up spectroscopy of their central stars with SALT. We have 
found that the spectra of both stars are very similar to 
those of the well-known LBVs P Cygni and AG Car (dur- 
ing its brightness minima in 1985-1990 and 2003) and the 
recently discovered cLBV MN112. We have also found that 
both stars, which we call WS1 and WS2, show significant 
photometric variability. Namely, WS1 brightened in the R- 
and /-bands by 0.68 ± 0.10 mag and 0.61 ± 0.04 mag, re- 
spectively, during the last 13-18 years, while WS2 (known 
as Hen 3-1383) varies its B,V,R,I and K B brightnesses by 
~ 0.5 — 0.9 mag on time-scales from 10 days to decades. 



9 Note that the fast-rotating runaway stars might also be the pro- 
genitors of long gamma-ray bursts (Hammer et al. 2006; Canticllo 
et al. 2007; Eldridge, Langer & Tout 2011; Le Floc'h et al. 2011). 



© 0000 RAS, MNRAS 000, 000-000 



14 V. V. Gvaramadze et al. 



Taken together, these findings strongly suggest that the de- 
tected objects are LBVs. The LBV classification of WS2 
is also supported by the spectral variability of this star on 
a timescale of ~ 3 months. To corroborate that WS1 and 
WS2 are bona fide LBVs, it is necessary to detect the ma- 
jor changes in their brightness (> 1 — 2 mag) and spectra. 
Further spectroscopic and photometric monitoring of these 
stars is therefore warranted. We have discussed a connection 
between the location of massive stars in the field and their 
fast rotation, and suggested that the LBV activity of WS1 
and WS2 might be directly related to their possible runaway 
status. 

To conclude, careful examination of the archival data 
of the WISE survey and the Spitzer Legacy programs will 
undoubtedly lead to new interesting discoveries, which hope- 
fully will advance our understanding the LBV phenomenon. 
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